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Abstract 

The physical nature of compressible turbulence is of fundamental importance in a variety of 
astrophysical settings. We present the first direct evidence that mean kinetic energy cascades con- 
\mJ servatively beyond a transitional "conversion" scale-range despite not being an invariant of the 

(-H compressible flow dynamics. We use high-resolution three-dimensional simulations of compress- 

1-H ible hydrodynamic turbulence on 512^ and 1024^ grids. We probe regimes of forced steady-state 

(-H isothermal flows and of unforced decaying ideal gas flows. The key quantity we measure is pressure 

^^ dilatation cospectrum, E^^{k), where we provide the first numerical evidence that it decays at a 

rate faster than k^^ as a function of wavenumber. This is sufficient to imply that mean pressure 
dilatation acts primarily at large-scales and that kinetic and internal energy budgets statistically 
decouple beyond a transitional scale-range. Our results suggest that an extension of Kolmogorov's 
inertial-range theory to compressible turbulence is possible. 
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1 Introduction 

Turbulence plays a critical role in essentially all astrophysical systems that involve gas dynamics. For 
example, interstellar turbulence is believed to be driven on large scales by differential galactic rotation 
or supernovae explosions and dissipated at the smallest scales by microphysical processes. It is widely 
believed that energy is transferred from the largest scales down to the dissipation scales through a cascade 
process. Measurements of interstellar scintillation (ISS) (1 ) caused by scattering in the interstellar 
medium shows a power-law electron density spectrum with a slope close to Kolmogorov's —5/3 scaling 
over 5 decades in scale. A composite spectrum combining ISS data with that of differential Faraday 
rotation angle and gradients in the average electron density also yields a Kolmogorov power-law over at 
least ten decades in scale, dubbed as "the great power-law in the sky" (|2])- High-resolution simulations 
of the interstellar medium ( 3 ) also reveal power-law scaling of spectra and structure functions which 
are reminiscent of incompressible turbulence (albeit with different slopes). Yet, observations of the 
interstellar medium suggest that the Mach number of turbulent motions is of order 0.1 to 10, such that the 
flow is often compressible ([!]). Similarly, high-resolution simulations of the intracluster medium (ICM) 
between galaxies ([5] IB]) reveal Kolmogorov-like spectra of kinetic energy at subsonic and transonic Mach 
numbers. It is not clear how Kolmogorov's 1941 phenomenology, which forms the cornerstone for our 
understanding of incompressible hydrodynamic turbulence, can carry over to such flows which exhibit 
significant compressibility effects. Several recent studies have addressed the problem of compressible 
hydrodynamic turbulence numerically (e.g. [71 [31 HI ISl HO]) and theoretically (e.g. [TTl [TH [T3] ) , but a 
physical description equivalent to that of Kolmogorov's still eludes us. 

The idea of a cascade itself is without physical basis since kinetic energy is not a global invariant of the 
inviscid dynamics. Compressible flows allow for an exchange between kinetic and internal energy through 
two mechanisms: viscous dissipation and pressure dilatation. While the former process is localized to 
the smallest scales just like in incompressible turbulence, the latter is a hallmark of compressibility and 
can a priori allow for an exchange at any scale through compression and rarefaction. Recently, [141 115] 
showed how the assumption that pressure dilatation cospectrum decays fast enough rigorously implies 
that mean kinetic energy cascades conservatively despite not being an invariant. The assumption entails 
that the mechanism of pressure dilatation acts primarily at the largest scales and vanishes on average at 
smaller scales. In this Letter, we shall provide the first empirical evidence in support of this assumption. 

The outline of this Letter is as follows. Section [2] describes our numerical simulations and Section 
[3] presents our coarse-graining method for analyzing nonlinear scale interactions. Section [4] discusses 
the significance of pressure dilatation cospectrum and the assumption we are testing. Our main results 
are then described in Section [5] followed by a discussion on the role of shocks in Section [6[ The Letter 
concludes with Section [71 



2 Numerical simulations 

We analyze data from four numerical simulations of compressible turbulence, summarized in Table [ll by 
solving the continuity, momentum, and internal energy equations 

dtp +V.(pu) = , (1) 

dt[pMi) +V-{puu) = -VP + pF , (2) 

dtipe) +V-(/9eu) = -PV-u , (3) 

supplemented with an equation of state (EOS) for the fluid. Here, u is velocity, p is density, P is 
pressure, e — P/[(7 — l)p] is internal energy per unit mass for a heat capacity ratio 7, and F is an 
external acceleration field stirring the fluid. We carry out two kinds of simulations; runs I and III in 
which the fluid is isothermal and the flow is constantly driven by a non-zero F, and runs II and IV in 
which the fluid is an ideal gas and the turbulence is decaying with F = (see Table fl]). 

The simulation domain is a periodic box T — [0, 27r)'^. For isothermal forced runs I and III, we start 
with a uniform density field p(x,t = 0) = I. The forcing function is adapted from 16 and has the form 
F(x, t) — J2i=i ^O^j) s^P (ikj'x) + complex conjugate. The Fourier amplitude, F(kj), has a compressive 
component parallel to kj and another solenoidal component perpendicular to kj. On average, the two 
components are of equal magnitude. The amplitudes and phases of F(kj) are random in time. The 
22 forced wavevectors, fcj, are nearly isotropically distributed in the range 2 < \kj\ < a/12. The initial 
conditions for unforced decaying Runs II and IV are taken from the last recorded output of Runs I and 
III, respectively, after reaching a statistically steady state. 

We use the central finite- volume scheme on overlapping cells ([17]) to solve eqs. (|l])-(l2]) in conservative 
form. Instead of solving internal energy eq. (p|, the code solves the equivalent total energy equation. 
The code can also solve compressible ideal magnetohydrodynamics. Details of the whole algorithm 
have been documented in [TH] and [H]. Our method has been verified to achieve the expected order of 
accuracy and to have very low numerical dissipation. The high-order, low-dissipation, and divergence- 
free properties of this method make it an excellent tool for simulating compressible hydrodynamic and 
magnetohydrodynamic turbulence. 

Figure [l] plots spectra of the solenoidal and compressive components of the velocity field from runs 
III and IV, both of which exhibit putative Kolmogorov-like power-law scaling over wavenumber range 
k e [2,50]. Here, the velocity field is decomposed as u(x) — u*(x) -I- u'^(x), such that u* is solenoidal 
and u"^ is irrotational. Spectra from run III are averaged over times after kinetic energy has reached 
steady state. Spectra from the decaying case, run IV, are normalized at each time by the value of (|up) 
at that time. Here, (...) is a volume average over the domain T, J^. d^x{. . .). 



Run 


N'-^ 


Flow 


EOS 


Mt 


I 


512^ 


forced 


isothermal 


0.43 


II 


5123 


decaying 


ideal gas 


- 


III 


10243 


forced 


isothermal 


0.44 


IV 


10243 


decaying 


ideal gas 


- 



Table 1: Simulation parameters: N^ is the grid-size, Mt = Urms/cth is the Mach number and ctn is the 
sound speed. In the decaying runs, initial Mt is similar to that of forced runs but monotonically decays 
in time. 





Figure 1: Plots of velocity spectra from run III (top) and run IV (bottom). Solid (blue) plot shows 
spectrum of the solenoidal velocity component. Dashed (red) shows compressive velocity component. 
Straight lines are for reference and have a — 5/3-slope. 

3 Analyzing nonlinear scale interactions 



The key analysis method we use is a "coarse-graining" or "filtering" approach common in Large Eddy 
Simulation (LES) literature on turbulence modelling. See for example [2011111 [22 and references therein. 
The approach was developed by [23l [24| to analyze nonlinear scale interactions in flow fields. It was 
further refined and utilized by [25j and extended to magnetohydrodynamic ([17]), geophysical (;2B]), and 
compressible ( [HI UHl [H] ) flows. The method itself is simple. For any field a(x), a "coarse-grained" or 
(low-pass) filtered field, which contains modes at length-scales > i, is defined as 



a^(x) = / d3r G'^(r)a(x + r), 



(4) 



where G(r) is a normalized convolution kernel, / cPr G{r) = 1. An example of such a kernel is the 
Gaussian function, G{r) ~ -j^e^"^ /^. Its dilation G^(r) = l^'^Giy: jl) has its main support in a ball of 
radius i. Operation Q may be interpreted as a local space average. 

From the dynamical equation of field a(x), coarse-grained equations can then be written to describe 
the evolution of a£(x) at every point x in space and at any instant of time. Furthermore, the coarse- 
grained equations describe flow at scales > ^, for arbitrary i. The approach, therefore, allows for the 
simultaneous resolution of dynamics both in scale and in space, similar to wavelet analysis, and admits 
intuitive physical interpretation of various terms in the coarse-grained balance. 

Moreover, coarse-grained equations describe the large-scales whose dynamics is coupled to the small- 
scales through so-called subscale or subgrid terms. These terms depend inherently on the unresolved 
dynamics which has been filtered out. The approach thus quantifies the coupling between different scales 
and may be used to extract certain scale-invariant features in the dynamics. 

3.1 Analyzing Compressible Flows 



[14l 128] showed how a Favre (or density- weighted) decomposition can be employed to extend the coarse- 
graining approach to compressible turbulence. A Favre filtered field is weighted by density as 

In the rest of this Letter, we shall take liberty of dropping subscript £ whenever there is no risk of 
ambiguity. The resultant large-scale dynamics for continuity and momentum are, respectively, 

dtp + di{pui) = 0. (6) 

dtput + dj {pui Uj) = - dj {p T{ui ,Uj)) ~ diP + pFi , 

(7) 

where 

pT{Ui, Uj) = p{lMUj - Ui Ui) (8) 

is the subgrid stress from the eliminated scales < £. It is also straightforward to derive a kinetic energy 
budget for the large-scales, which reads 

dtP^^ + V-3i = -n, - A, + PiV-Ui + e"^' , (9) 

where Jf(x) is spatial transport flux of large-scale kinetic energy, and e™''(x) is the energy injected due 
to external stirring. The definition of these terms can be found in [T31[Tn]. — PV-u is large-scale pressure 



dilatation, and H^ (x) + A^(x) is the subgrid scale (SGS) kinetic energy flux to scales < I, 

nf(x) = -pdjUiT{ui,Uj) (10) 

A,(x) = -d,PT{p,u,) (11) 

where 

rip,u) = {pu)f~Pi,ue (12) 

in expression (11) is the subgrid mass- flux. Equations ([6|-(|9]) describe the dynamics at scales > £, for 
arbitrary i, at every point x and at every instant in time. They hold for each realization of the flow 
without any statistical averaging. 

The SGS flux is comprised of deformation work, Hi, and baropycnal work, A^, which are discussed 
in some detail in |28j . These represent the only two processes capable of direct transfer of kinetic energy 
across scales. Pressure dilatation, — P^V-u^, does not contain any modes at scales < i (or a moderate 
multiple thereof), at least for a filter kernel G(k) compactly supported in Fourier space. Therefore, 
pressure dilatation cannot participate in the inter-scale transfer of kinetic energy and only contributes 
to conversion of large-scale kinetic energy into internal energy. This was a crucial observation made in 
[H] on which the analysis in this Letter will be based. 

4 Pressure dilatation cospectrum 

Kinetic and internal energy budgets couple through two mechanisms. The first is viscous dissipation 
which was proved in [55] to be confined to the smallest scales i < i^ (the dissipation scale-range). 
Therefore, large-scale kinetic energy in Q docs not couple to internal energy via viscous dynamics for 
i 3> ip,. The second mechanism is pressure dilatation, — PV-u, which exchanges kinetic and internal 
energy via compression and rarefaction. 

It was shown in [14_ that if pressure dilatation co-spectrum, defined as 

E^^ik)^ ^ -P(k)V:u(-k), (13) 

fc-0.5<|k|<A;-|-0.5 

decays fast enough as a function of wavenumber, 

\E''''{k)\ < Cu™, Prras {kL)-^, /? > 1, (14) 

then mean pressure dilatation exchanges mean kinetic and internal energy over a transitional "conver- 
sion" scale-range of limited extent. Here, C is a dimensionless constant and L is an integral scale. At 
smaller scales beyond the conversion range, mean kinetic and internal energy budgets statistically decou- 
ple giving rise to an inertial range over which mean kinetic energy undergoes a scale-local conservative 




Figure 2: Magnitude of pressure dilatation cospectrum from run III (top) and run IV (bottom). Straight 
lines are for reference and extend over the fitting range. Top: straight lines have slopes of —1 and —2; 
power-law fit of data gives a slope of -1.88. Bottom: straight lines have slopes of —1 and —5/3; power-law 
fit of data gives a slope of -1.61. 



cascade. In this Letter, we provide the first empirical evidence in support of assumption (14) in Figure 
[2] which we shall discuss more below. 



The idea behind assumption ( 14 ) is straightforward and rests on the convergence of a series or an 



integral at infinity. In the limit of large Reynolds number, assumption ( 14 1 implies that mean large-scale 
pressure dilatation, PD{i) = — (P^V-Uf), asymptotes to a finite constant, 6 = — (PV-u), as £ ^- 0. In 
other words, PD{i) acting at scales > £ converges and becomes independent oi £ at small enough scales: 



lim PDie)= lim V E^^(k)::^ 



(15) 



^<k<K 



Note that PD(£) in (151 (also shown in Figure^ is a cumulative quantity representing the contribution 
from all wavenumbers k < K ~ i~^, whereas spectra E^^{k) and E^ik) in Figures 111 ^ are density 



functions. Convergence of PD{£) in (15) expresses the decoupling (in an average sense) between large- 
scale kinetic and internal balances. Such a decoupling is statistical and does not imply that small scales 
evolve according to incompressible dynamics. However, while small-scale compression and rarefaction 
can still take place pointwise, they yield a vanishing contribution to the space-average. 
We denote the smallest wavenumber at which such statistical decoupling occurs by Kc- 
Over the ensuing wavenumber-range, Kc < K <t^ K^ (where K^ = i^^), net pressure dilatation 
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Figure 3: Solid (blue) plot shows mean SGS kinetic energy flux, (11;^: + t^K)-, and dashed (red) plot 
shows mean pressure dilatation, PD{K), from run III (top) and run IV (bottom). 



does not play a role. For £ = K ^ over this wavenumber-range, budget tel) of large-scale kinetic energy 
becomes 



(n, 



Af) = 



(16) 



in a steady state and after space-averaging. If (e™-') is localized to the largest scales as discussed in [28], 
then {Hg + A() will be a constant, independent of scale £. 

A constant SGS flux implies that mean kinetic energy cascades conservatively to smaller scales, 
despite not being an invariant of the governing dynamics. This was one of the main conclusions in 
[14j . In particular, kinetic energy can only reach dissipation scales via the SGS flux, II^ + A^, through 
a scale-local cascade process. We are therefore justified in calling wavenumber-range Kc < K <^ K^ 
the inertial range of compressible turbulence. The existence of this inertial range in such flows warrants 
expectations that spectra with power-law scalings should exist and that their observation is evidence of 
a turbulent cascade process similar (although not necessarily identical) to that in incompressible flows. 

5 Results 

We will now present results from our simulations to provide the first empirical evidence in support of 
assumption (14|. Figure^shows that \E^^{k)\ decays at a rate significantly faster than fc^-'^, well in 
agreement with condition (|14|). We fitted the measurements with power laws over wavenumber range 



k £ [2,50]. We obtain a scaling of \EP'^{k)\ - k'^-^^ from run III and \EP'^{k)\ - k''^-^^ from run 
IV. The measurements shown from runs III and IV are averaged over 21 and 28 evenly-spaced time 
snapshots, respectively. The plots of \E^^{k)\ from run IV are normalized at each time by the value of 
^j, \E^-^{k)\ at that time. We have verified that results from our 512^ runs I and II are consistent with 
those presented. 

To further illustrate the main idea of this Letter, we analyzed the SGS flux terms, (n^+A^), in budget 
([9]). Figure [3] shows that pressure dilatation, PD{£), tends to a constant beyond a transitional conversion 
range. Since PD{£) is a cumulative quantity representing the contribution from all wavenumbers k < 
K — £~^ , its convergence to a constant beyond K — 10 implies that any smaller scales with fc > 10 give a 
negligible contribution to pressure dilatation. The SGS kinetic energy flux also becomes approximately 
const an1r] beyond this conversion range, implying that kinetic energy is cascading conservatively (without 
any "leakage" to/from internal energy) until the dissipation range is reached. It is notoriously hard to 
achieve a constant SGS flux range in limited resolution simulations due to viscous contamination. Even 
the largest-to-date simulation of incompressible turbulence on a 4096'^ grid by [3D] does not exhibit a 
clear constant flux. However, the plots in Figure [3] suggest that an inertial range over which SGS flux is 
constant does indeed arise at smaller scales over which PD{£) plateaus. 

6 The role of shocks 

Our results and the physical picture we are advancing might seem counter-intuitive at first. After all, a 
hallmark of compressible turbulence is the formation of shocks and the generation of sound waves. Such 
phenomena involve compression and rarefaction at all scales and are not restricted to small wavenumbers. 
However, there is no contradiction between the existence of such phenomena and our conclusions. 
Our results concern global pressure dilatation, — (PV-u), and not the pointwise quantity. While we 
expect very large pressure dilatation values in the vicinity of small-scale shocks, our results indicate 
that such a contribution will vanish when averaging over the flow domain due to cancellations between 
compression and rarefaction regions. For example, linear (small-amplitude) sound waves yield zero 
mean pressure dilatation, — (PV-u) = 0. Figure [4] shows that values of pointwise pressure dilatation 
at small scales, ~[PV-u(x) — Pi'V-Uf{x)], are an order of magnitude more intense relative to large- 
scale pressure dilatation, — P£V'Uj»(x), as expected. The distribution of — [PV'u(x) — Pi'V-ue{x.)] has 
heavy tails implying spatially rare but intense two-way exchange between kinetic and internal energy. 
Furthermore, the small-scale distribution has positive skewness implying that conversion from kinetic 
into internal energy through compression occupies less volume and is more intense relative to conversion 



^Unlike the cumulative quantity PD{£) which represents the contribution from all wavenumbers k < K = £~^, the 
SGS flux (n^ -I- Ai) is a hybrid quantity involving both wavenumbers k < K = i~^ (large-scale terms) and k > K = £~^ 
(subgrid terms) as definitions (|10|,Jll|l show. 



from internal into kinetic energy through rarefaction. Despite the larger intensity of conversion at small 
scales, its global contribution, — (PV-u — P^V-u^) = 2.0 x 10~^, is negligible compared to the mean of 
large-scale conversion, PD{i) = — (P^V-u^) = 3.7 x 10^^, for (r^ = K = 20. This is consistent with 
the plateau of PD{i) in Figure [s] 

The origin of such cancellations at small-scales may be understood through the following argument 
(|14j). While pressure in — (PV-u) derives most of its contribution from the largest scales, V-u is 
dominated by the smallest scales in the flow. Therefore, pressure varies slowly in space, primarily at 
scales ~ L, while V-u varies much more rapidly, primarily at scales £^ <C i, leading to a decorrelation 
between the two factors. Our numerical results presented in this Letter can be regarded as empirical 
support to such a physical argument. 

7 Conclusions &; Discussion 

The main result of our Letter was to provide the first empirical evidence that pressure dilatation co- 



spectrum decays at a rate faster than fc~^, in accord with condition (14 1. This is sufficient to imply that 
exchange of kinetic and internal energy through compression and rarefaction takes place at the large 
scales, on average. At smaller scales beyond a transitional "conversion" scale-range, mean kinetic and 
internal energy budgets statistically decouple and kinetic energy can only reach dissipation scales via a 
scale-local conservative cascade process. The existence of such a conservative cascade in compressible 
turbulence justifies expectations that spectra with power-law scalings should exist in such flows. Further- 
more, it suggests that observations of power-law spectra in astrophysical systems are indeed evidence of 
a turbulent cascade process similar (although not necessarily identical) to that in incompressible flows. 
This work leads us to some new questions which we hope to address in future studies. Does the 
decay rate of pressure dilation cospectrum depend on the decaying/forced nature of the turbulence or on 
other parameters of the flow such as equation of state, Mach number, and the ratio of compressive-to- 
solenoidal components of the velocity field? Is the power law scaling of pressure dilation cospectrum that 
we find in Figure [2] reflecting the asymptotic scaling at arbitrarily high Reynolds numbers? Does the 
value of PD{£) in Figure [3] relative to that of SGS flux increase by increasing the compressive amplitude 
of the forcing? We invite future studies to investigate these questions and to try and reproduce our 



results presented in this Letter. Verifying assumption (14) under a variety of controlled conditions 
would substantiate the idea of statistical decoupling between kinetic and internal energy budgets. This 
is of paramount importance for future attempts to extend Kolmogorov's ideas on a conservative cascade 
of kinetic energy to compressible turbulence. 

Acknowledgements. We are grateful to J. Cho for providing us with his forcing subroutine. H. A. thanks 
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Figure 4: For i^^ ^ K ~ 20 in Fig. ^from run IV, top two panels visualize pointwise pressure dilatation 
at large scales (top-left), —Pi'V-ue{x.), and the residual from small-scales (top-right), — PV'u(x) -|- 
PfV'Uf(x), in a 1024^ domain cross-section. The color maps are rescaled to fit within the color bars' 
bounds. Bottom two panels show the distributions of large-scale (bottom-left) and small-scale (bottom- 
right) pressure dilatation. Dashed curves are for reference and show unnormalized gaussian distributions 
with zero mean. 
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